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Using a recently developed procedure - multiple wave packet decomposition - here we study 
the phase time formulation for tunneling/reflecting particles colliding with a potential barrier. 
To partially overcome the analytical difficulties which frequently arise when the stationary phase 



' method is employed for deriving phase (tunneling) time expressions, we present a theoretical exer- 

o ' 

cise involving a symmetrical collision between two identical wave packets and an one-dimensional 



rectangular potential barrier. Summing the amplitudes of the reflected and transmitted waves 
- using a method we call multiple peak decomposition - is shown to allow reconstruction of the 



scattered wave packets in a way which allows the stationary phase principle to be recovered. 
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Recently, a series of experimental results 



flflflQ, 



some of them confirming the pos- 



sibility of superluminal tunneling speeds for photons, have revived an interest in the tun- 
neling time analysis , |9j. On the theoretical front, people have tried to introduce 
quantities that have the dimension of time and can somehow be associated with the pas- 
sage of the particle through the barrier or, strictly speaking, with the definition of the 



;unneling^ time. Th e se proposals 
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13, 



rave led to the introduction of several time definitions 

3, Q, 



some of which can be organized into three 
groups. (1) The first group comprises a time- dependent description in terms of wave pack- 
ets where some features of an incident packet and the comparable features of the transmitted 
packet are utilized to describe a quantifiable delay as a tunneling time 9j. (2) In the sec- 
ond group the tunneling times are computed based on averages over a set of kinematical 
paths, whose distribution is supposed to describe the particle motion inside a barrier. In 
this case, Feynman paths are used like real paths to calculate an average tunneling time 
with the weighting function exp [i S x(t)/h], where S is the action associated with the path 
x(t) (where x(t) represents the Feynman paths initiated from a point on the left of the 
barrier and ending at another point on the right of it [2l|]). The Wigner distribution paths 



3, 



and the Bohm approach 



22 



231 ] are included in this group. (3) In the third group 



we notice the introduction of a new degree of freedom, constituting a physical clock for the 
measurements of tunneling times. This group comprises the methods with a Larmor clock 
ll| or an oscillating barrier 24J. Separately, standing on itself is the dwell time defined 
by the interval during which the incident flux has to exist and act, to provide the expected 
accumulated particle storage, inside^the barrier 



There is no general agreement [5|, [8( among the above definitions about the meaning of 
tunneling times (some of the proposed tunneling times are actually traversal times, while 
others seem to represent in reality only the spread of their distributions) and about which, 
if any, of them is the proper tunneling time 5] . In the context where we intend to work on, 
the tunneling mechanism is embedded by theoretical constructions involving analytically- 
continuous gaussian, or infinite-bandwidth step pulses to examine the tunneling process. 
Nevertheless, such holomorphic functions do not have a well-defined front in a manner that 
the interpretation of the wave packet speed of propagation becomes ambiguous. Moreover, 
infinite-bandwidth signals cannot propagate through any real physical medium (whose trans- 
fer function is therefore finite) without pulse distortion, which also leads to ambiguities in 
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determining the propagation velocity during the tunneling process. For instance, some of 
the barrier traversal time definitions lead, under tunneling time conditions, to very short 
times, which can even become negative. It can precipitately induces an interpretation of 
violation of simple concepts of causality. Otherwise, negative speeds do not seem to create 
problems with causality, since they were predicted both within special relativity and within 
quantum mechanics 18J. A possible explanation of the time advancements related to the 
negative speeds can come, in any case, from consideration of the very rapid spreading of the 
initial and transmitted wave packets for large momentum distribution widths. Due to the 
similarities between tunneling (quantum) packets and evanescent (classical) waves, exactly 
the same phenomena are to be expected in the case of classical barriers [47]]. The existence 
of such negative times is predicted by relativity itself based on its ordinary postulates , 



and they appear to have been experimentally detected in many works 
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23]. 



In this extensively explored scenario, the first group quoted above contains the so-called 



phase times 
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29 



30l | which are obtained when the stationary phase method (SPM) 3^] is 



employed for obtaining the times related to the motion of the wave packet spatial centroid. 
Generically speaking, the SPM essentially enables us to parameterize some subtleties of 



onances 
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as well as the Hartman effect |3 



several quantum phenomena, such as tunneling 2j, |9 L 
reflection and incidence-transmission interferences 
and its superluminal traversal time interpretation [5|, ]7J, 
most usual approximation method for describing the group velocity of a wave packet in a 
quantum scattering process represented by the collision of a particle with a potential barrier 



19]. In fact, it is the simplest and 
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39]. 



In the following study we will concentrate on some incompatibilities that appear when the 
SPM is utilized for deriving tunneling times. After quantifying the restrictive conditions for 
the use of the method, at the end of our analysis, we discuss a theoretical exercise involving 
a symmetrical collision between two identical wave packets and an one- dimensional rectan- 
gular potential barrier. We demonstrate that by summing the amplitudes of the reflected 
and transmitted waves in the scope of what we denominate a multiple peak decomposition 
analysis [39]], we can recompose the scattered wave packets in a way that the analytical 
conditions for the SPM applicability are totally recovered. 

The SPM can be successfully applied for describing the movement of the center of a wave 
packet constructed in terms of a momentum distribution g(k — k ) which has a pronounced 



3 



peak around k . By assuming that the phase that characterizes the propagation varies 
smoothly around the maximum of g(k — k ), the stationary phase condition enables us to 
calculate the position of the peak of the wave packet (highest probability region to find the 
propagating particle). With regard to the tunneling phenomenon, the method is usually 
applied to find the position of a wave packet that traverses a potential barrier. For the case 
in which we consider a rectangular potential barrier V(x), V(x) — V if x £ [-L/2, L/2] 
and V(x) = if x 4 [-L/2, L/2], 

V x € [-L/2, L/2] 
VU) = \ (1) 
x $ [-L/2, L/2] 

it is well known that the transmitted wave packet solution (x > L/2) calculated by means 



(2) 



of the Schroedinger formalism is given by 
f w dk 

ip T (x,t) = J — g(k-k )\T(k,L)\exp 

In case of tunneling, the transmitted amplitude and the phase shift are respectively given 
by 



ikix- L/2) -i t + i @(k, L) 

2 m 



and 



\T(k,L)\ 



S(k,L) 



1 + 



w 



arctan 



4 k 2 p 2 (k) 
2 k 2 



sinh 2 [p(k) L] 
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kp(k) 



tanh [p(k) L] 
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(3) 



(4) 



for which we have made explicit the dependence on the barrier length L, and we have 



adopted p{k) 



[w 



k 2 ) 2 with w = (2mV ) 2 and h = 1. By not considering any eventual 



distortion that \T(k,L)\ could cause to the supposedly symmetric function g{k — k ), the 
stationary phase condition is indiscriminately applied to the phase leading to 
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dk 



k 2 

k(x-L/2) t + Q(k,L 



2 m 







x _ L /2-~t + 



dQ{k,L) 



m 



dk 



0. 



(5) 



The above result is frequently adopted for calculating the transit time tx of a transmitted 
wave packet when its peak emerges at x — L/2, 



tr 



m 



dQ(k,a {L) ) 



k„ 



dk 



2 m L j sinh (a) cosh (a) — (2 k 2 ^. 



w 2 ) k 2 ^ a 



'"max 



4fcLx {w 2 -klJ+w^ sinh 2 (a) 



(6) 



where we have defined the parameter a = (w — k^J) 2 L. The concept of opaque limit is 
introduced when we assume that fe max is independent of L and then we make a tend to oo 

[3]. In this case, the transit time can be rewritten as 

ol 2 m , . 

tT = k p(k V (7) 

'•'max /-*V max / 

In the literature, the value of fc max is frequently approximated by k , the maximum of g(k—k ), 
which, in fact, does not depend on L and could lead us to the superluminal transmission 



time interpretation 
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461 ] . To clear up this point, we notice that when we take the so 
called opaque limit in Eq. ([7j) , with L going to oo and w fixed as well as with w going to oo 
and L fixed, with k < w in both cases, the expression (J7j) leads to times corresponding to 



a transmission process performed with velocities larger than c 



19]. 



Such a superluminal interpretation was extended to the study of quantum tunneling 
through two successive barriers separated by a free region j^]. In this approach, the total 
traversal time should be independent of the barrier widths and of the distance between 
the barriers. In a subsequent analysis, the same technique was applied to a problem with 
multiple successive barriers where the tunneling process was designated as a highly non-local 



phenomenon 



46|]. 



It would be perfectly acceptable to consider fc max = k for the application of the stationary 
phase condition if the momentum distribution g(k — k ) centered at k was not modified by 
any boundary condition. That is the case of the incident wave packet before the collision 
with the potential barrier. In this sense, and in the context of the above quoted theoretical 
results, our criticism is concerned with the way of obtaining all the above results for the 
transmitted wave packet. It has not taken into account the bounds and enhancements 
imposed by the analytical form of the transmission coefficient. 

To perform the correct analysis, we should calculate the correct value of fc max to be 
substituted in Eq. ([6]) before taking the opaque limit. We are thus obliged to consider the 
relevant amplitude for the transmitted wave as the product of a symmetric momentum 
distribution g(k — k ), which describes the incoming wave packet, by the modulus of the 
transmission amplitude T(k,L), which is a crescent function of k. The maximum of this 
product representing the transmission modulating function would be given by the solution 
of the equation 



g(k-k ) \T(k,L) 



g'(k-k ) + \T(k,L) 



g(k-k ) \T(k,L)\\ 



0. (8) 



Obviously, the peak of the modified momentum distribution is shifted to the right of k 
so that /c max has to be found in the interval ]/c ,w[. Moreover, we can demonstrate by the 
numerical results of Table 1 that fc max presents an implicit dependence on L. For obtaining 
the Table 1 data we have found the maximum of g(k — k ) \T(k, L) \ by assuming a gaussian 
distribution g(k — k () ) = (f^J ex P — a ( k ~ k °) almost completely comprised in the interval 
[0,w\. 

By increasing the value of L with respect to the wave packet width a, the value of fc max 
obtained from the numerical calculations to be substituted in Eq. ([6]) also increases up to L 
reaches certain values for which the modified momentum distribution becomes unavoidably 
distorted. In this case, the relevant values for k are concentrated in the neighborhood of the 
upper boundary value w. We shall show in the following that the value of L, which sets up 
the distortion the momentum distribution can be analytically obtained in terms of a. 

Now, if we take the opaque limit of a by fixing L and increasing w, the above results 
immediately ruin the superluminal interpretation upon the result of Eq. (jH]), since t^ L tends 
to oo when fc max is substituted by w. Otherwise, when w is fixed and L tends to oo, the 
parameter a calculated at k = w becomes indeterminate. The transit time still tends to 
oo but now it exhibits a peculiar dependence on L, which can be easily observed by defining 
the auxiliary function 

sinh (a) cosh (a) — a 

G{a) = r-TTT-s • (9) 

sinn [a) 

When «>1, the transmission time assumes infinite values 

2mL 2m 
trp = G(a) =>• t T ~ j > oo. (ID) 

with an asymptotic dependence on (w 2 — k 2 ) Only when a tends to we have an explicit 
linear dependence on L given by 

2mL v _ \ G(a) \ 4mL 



t% = lim <^ } = (11) 

w »->o { a J 3w 

In addition to the above results, the transmitted wave must be carefully studied in terms 
of the ratio between the barrier extension L and the wave packet width a. For very thin 
barriers, i. e. when L is much smaller than a, the modified transmitted wave packet presents 
substantially the same form of the incident one. For thicker barriers, but yet with L < a, 
the peak of the gaussian wave packet modulated by the transmission coefficient is shifted 
to higher energy values, i. e. fc max > k increases with L. For very thick barriers, i. e. 



TABLE I: The values of k numerically obtained in correspondence with the increas- 
ing of the barrier extension L. The values are calculated in terms of the wave packet 
width a for different values of the potential barrier height expressed in terms of w a. 
We have fixed the incoming momentum by setting k a = 1. 



w a 

L/a 




1.5 




2.0 




4.0 




6.0 




8.0 




10 




20 


0.00 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


0.10 


1 


0235 


1 


0648 


1 


3799 


1 


6769 


1 


8547 


1 


9397 


2 


0051 


0.20 


1 


0794 


1 


1825 


1 


6571 


1 


9178 


2 


0000 


2 


0204 


2 


0203 


0.30 


1 


1478 


1 


3001 


1 


8430 


2 


0289 


2 


0562 


2 


0551 


2 


0342 


0.40 


1 


2196 


1 


4116 


1 


9874 


2 


1025 


2 


0986 


2 


0857 


2 


0484 


0.50 


1 


2921 


1 


5194 


2 


1155 


2 


1668 


2 


1399 


2 


1170 


2 


0628 


0.60 


1 


3649 


1 


6266 


2 


2429 


2 


2314 


2 


1828 


2 


1495 


2 


0775 


0.70 


1 


4383 


1 


7360 


2 


3819 


2 


3002 


2 


2281 


2 


1834 


2 


0925 


0.80 




* a 


1 


8489 


2 


5466 


2 


3751 


2 


2761 


2 


2188 


2 


1078 


0.90 




* 


1 


9646 


2 


7627 


2 


4578 


2 


3272 


2 


2558 


2 


1234 


1.00 




* 




* 


3 


1137 


2 


5504 


2 


3818 


2 


2947 


2 


1392 



Tor the values of L marked with *, we can demonstrate by means of Eqs. (|12M13[) that the 
modulated momentum distribution has already been completely distorted. In this case, the 
maximum has no meaning in the context of the applicability of the method of stationary 
phase. 

when L > a, we are able to observe that the form of the transmitted wave packet is badly 
distorted with the greatest contribution coming from the Fourier components corresponding 
to the energy w just above the top of the barrier in a kind of filter effect. We observe that 
the quoted distortion starts to appear when the modulated momentum distribution presents 
a local maximal point at k = w which occurs when \g{k — k ) \T(k,L)\]\,_ > 0. Since 
the derivative of the gaussian function g{k — k ) is negative at k = w, the previous relation 
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gives 

~T'(k,L) 



g'iw - fc ) 

— ; z-t- < lim 

g(w — K ) k^w 



T{k, L) 



1 + hJ±£ 
1 ^ 4 



< (12) 



which effectively represents the inequality 



T^-^X— L> V2 a ( -J' (3) 

Due to the /z/£er effect, the amplitude of the transmitted wave is essentially composed 
by the plane wave components of the front tail of the incoming wave packet that reaches 
the first barrier interface before the peak arrival. Meanwhile, only whether we had cut the 
momentum distribution off at a value of k smaller than w, i. e. k ~ (1— 5)w, the superluminal 
interpretation of the transition time (J7J) could be recovered. In this case, independently of 
the way that a tends to oo, the value assumed by the transit time would be approximated 
by £ j. ~ 2m I w 5, which is a finite quantity. Such a finite value would confirm the hypothesis 
of superluminality. However, the cut off at k « (1 — S)w increases the amplitude of the 
tail of the incident wave as we can observe in Fig. HJ It means that the contribution of 
wave packet tail for the final composition of the transmitted wave is put on the same level 
with the contribution of the peak of the incident wave. Consequently, an ambiguity in the 
definition of the arrival time is created. 

To summarize, at this point we are particularly convinced that the use of a step- 
discontinuity to analyze signal transmissions in tunneling processes deserves a more careful 
analysis than the immediate application of the stationary phase method. The point is that 
we cannot find an analytic-continuation between the above-barrier case solutions and the 
below-barrier case solutions. By assuming the factual influence of the amplitude of the trans- 
mitted wave, we may introduce an alternative analysis where we consider the possibility of 
using the multiple peak decomposition technique developed for the above barrier diffusion 
problem [39(. By means of such an experimentally verifiable exercise, we shall be able to 
understand how the filter effect can analytically affect the calculations of transit times in 
the tunnel process. 

In the framework of the multiple peak decomposition 39|, we suggest a suitable way for 
comprehending the conservation of probabilities for a very particular scattering configura- 
tion where the asymmetric aspects above discussed can be totally eliminated. In order to 
recover the scattered momentum distribution symmetry conditions for accurately applying 
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FIG. 1: Dependence of the wave packet shape on the cut off value of a momentum distribution 
centered around k Q = 0.5w with the values of k comprised between and k cut off . 

the SPM, we assume a symmetrical colliding configuration of two wave packets traveling 
in opposite directions. By considering the same rectangular barrier V(x), we solve the 
Schroedinger equation for a plane wave component of momentum k for two identical wave 
packets symmetrically separated from the origin x — 0. At time t = —(mL) / (2k ) chosen for 
mathematical convenience, we assume that they perform a totally symmetric simultaneous 
collision with the potential barrier. The wave packet reaching the left (right) side of the 
barrier is represented by 

+ 00 

^ L ( R )(x,t)= dkg(k- k o )0 L{R \k,x) exp[-iEt]. (14) 
•'o 

Here we have assumed that the limits of the above integral can be naturally extended from 
the interval [0, w] to the interval [0, oo] as a first approximation. Its range of validity can 
be controlled by the choice of the width Ak of the momentum distribution g{k — k ) (with 
k > 0) with Ak enhanced by the barrier's height (V ). By assuming that <f) L( - R) {k,x) are 

9 



(f) L ( R Xk,x) 



Schroedinger equation solutions, at the time t = —(mL)/(2k ), i. e. when the wave packet 
peaks simultaneously reach the barrier, we can write 

p\ {R \k,x) = exp[±ikx] + R^ {R \k,L)exp[Tikx] x < -L/2 {x > L/2), 

pV R) (k,x) = aY R \k)exp[Tpx] + /% {R \k) exp [±p x] - L/2 < x < L/2, 
pY R \k,x) = T2 (R) (k,L)exp[±ikx} x > L/2 (x < -L/2). 

where the upper(lower) sign is related to the index L(R). By assuming the conditions for the 
continuity of <$^ ,R and their derivatives at x = —L/2 and x = L/2, after some mathematical 
manipulations, we can easily obtain 

and 

Tb R (k, L) = exp [_i k L] { exp[p(fc) L][l - exp[2 .e(fc,L)]] j 

where Q(k,L) is given by the Eq. (j3J) and B^{k,L) and T R {k,L) as well as R R {k,L) and 
T~{k, L) are intersecting each other. By analogy with the procedure of summing amplitudes 



that we have adopted in the multiple peak decomposition scattering [39(], such a pictorial 
configuration obliges us to sum the intersecting amplitude of probabilities before taking their 
squared modulus in order to obtain 

nL.ji/. n . r r R,L n T , r . , rl / exp[p(fc)L] +exp[iQ(k,L)} \ 

R B (k,L)+T B (k, L) = exp [~tkL]{ — — ^ £] — ^ ^ j = exp L + <p(k, L)]} 

(17) 



with 



V {k, L) = arctan < 2 kp(k) sinh [p(fc) L] 

' u> 2 + (& 2 - p 2 (£;)) cosh [p(k) L] 1 



From Eq. (IT7I) . it is important to observe that, differently from the previous standard tun- 
neling analysis, by adding the intersecting amplitudes at each side of the barrier, we keep the 
original momentum distribution undistorted since \B^' R (k, L) + T R ^(k, L)\ is equal to one. 
At this point we recover the most fundamental condition for the applicability of the SPM. 
It allows us to accurately find the position of the peak of the reconstructed wave packet 
composed by reflected and transmitted superposing components. The phase time interpre- 
tation can be, in this case, correctly quantified in terms of the analysis of the new phase 
ip(k,L). By applying the stationary phase condition to the recomposed wave packets, the 
maximal point of the scattered amplitudes g(k — k )\R^' R (k, L) +T R '^(k, L) \ are accurately 
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given by fc max = k so that the traversal/reflection time or, more generically, the scattering 
time, results in 



2m L w 2 sinh (a) — a k 2 

(19) 



k a 2 kl — w 2 + w 2 cosh (a) 

with a previously defined. It can be said metaphorically that the identical particles rep- 
resented by both incident wave packets spend a time of the order of tf T inside the barrier 
before retracing its steps or tunneling. In fact, we cannot differentiate the tunneling from 
the reflecting waves for such a scattering configuration. The point is that we have introduced 
the possibility of improving the efficiency of the SPM in calculating reflecting and tunneling 
phase times, by studying a process where the conditions for applying the method are totally 
recovered. We have demonstrated that the transmitted and reflected interfering amplitudes 
results in a unimodular function which just modifies the envelop function g(k — k ) by an 
additional phase. The previously pointed out incongruities which cause the distortion of the 
momentum distribution g(k — k () ) are completely eliminated in this case. At the same time, 
one could argue about the possibility of extending such a result to the tunneling process 
established in a standard way. We should assume that in the region inside the potential 
barrier, the reflecting and transmitting amplitudes should be summed before we compute 
the phase changes. Obviously, it would result in the same phase time expression as repre- 
sented by ( TT9|) . In this case, the assumption of there (not) existing interference between the 
momentum amplitudes of the reflected and transmitted waves at the discontinuity points 
x = —L/2 and x = L/2 is purely arbitrary. Consequently, it is important to reinforce the 
argument that such a possibility of interference leading to different phase time results is 
strictly related to the idea of using (or not) the multiple peak (de) composition in the region 
where the potential barrier is localized. To illustrate the difference between the standard 
tunneling phase time t? and the alternative scattering phase time t^ we introduce the new 
parameter n = k^Jw 2 and we define the classical traversal time r = (ml)/fc mai . Then we 
can obtain the rates 

tx 2 f cosh (a) sinh (a) — a n (2n — 1) 1 v tf T 2 f n a + sinh (a) 1 

Rt{oi) = - = - | [4n(1 _ n)+sinh2(a)] ) and «r(«) = - = - { 2n _ 1 + cosh(a) j 

(20) 

which are plotted in the Fig.([2]) for some discrete values of n varying from to 1., The most 
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FIG. 2: Time rates for the standard tunneling and the new scattering process. The rates R(a) 
and R^(a) can be understood as transmitted times in the units of the classical propagation time 
r. Both present the same asymptotic behavior, which, in a totally restrictive mathematical sense, 
and in the stationary phase analysis context, offers the possibility of a superluminal interpretation 
for the peak of the transmitted wave packet so that the SPM can be accurately applied. 

common limits of the above expressions are given by 

lim = km {R T (a)} = 0, (21) 

and 

lim {R T (a)} = 1 + — , and lim \Rl(a)\ = 1 + - (22) 
a^O 2n a^O L ' n 

Both present the same asymptotic behavior, which, in a totally restrictive mathematical 
sense of the stationary phase analysis context, allows the possibility of a superluminal inter- 
pretation for the peak of the transmitted wave packet. The main point is that, by now, from 

12 



the point of view of the analytical limitations, the SPM can be accurately applied. At this 
point it is convenient to notice that the superluminal phenomena, observed in the experi- 

4|, generated a 



0,0,3 



ments with tunneling photons and evanescent electromagnetic waves 
lot of discussions on relativistic causality. In fact, superluminal group velocities in connec- 
tion with quantum (and classical) tunnelings were predicted even on the basis of tunneling 
time definitions more general than the simple Wigner's phase time |29f] . Olkhovsky et al. 
discuss a simple way of understanding the problem J3|]. In a causal manner, it might consist 
in explaining the superluminal phenomena during tunneling as simply due to a reshaping 
of the pulse, with attenuation, as already attempted (at the classical limit) [36(. The later 
parts of an incoming pulse are preferentially attenuated, in such a way that the outcoming 
peak appears shifted towards earlier times even if it is nothing but a portion of the incident 
pulse's forward tail 2], In particular, we do not intend to expand on the delicate ques- 
tion whether superluminal group velocities can sometimes imply superluminal signaling. It 
is a controversial subject which has been extensively explored in the literature about the 
tunneling effect (jlj] and references therein). 

Turning back to the scattering time analysis, we can observe an analogy between our 
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The HE is 



results and the results interpreted from the Hartman Effect (HE) analysis 
related to the fact that, for opaque potential barriers, the mean tunneling time does not 
depend on the barrier width. For large barriers the effective tunneling velocity can become 
arbitrarily large so that the tunneling phase time becomes independent of the barrier width. 
It seems that the penetration time, needed to cross a portion of a barrier, in the case 
of a very long barrier starts to increase againafter the plateau corresponding to infinite 



speedproportionally to the distance [48j. Our phase time dependence on the barrier width 
is similar to that which leads to Hartman interpretation as we can infer from Eqs. (I21H22I) . 
Only when a tends to we have an explicit linear time dependence on L, 

,„ 2 m L ( 1 \ . . 

4 = 1 + - > 23 

w \ n ) 

which agrees with calculations based on the simple phase time analysis where = 
(l + if") ' as we can observe in the by Eq. (fTTj) for n = 1. 
Here at once it is important to emphasize that the wave packets for which we compute 
the phase times illustrated in the Fig.(T5]) are not constructed with the same momentum dis- 
tributions. The phase Q(k, L) appears when we treat separately the momentum amplitudes 
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g(k — k ) \T(k, L)\ and g(k — k )\R(k, L)\ and the other one (p(k, L) appears only when we 
sum the amplitudes g(k — k ) \T(k, L) + R(k, L)\ = g(k — k ) in order to obtain a symmetrical 
distribution. It requalifies the SPM for accurately computing the time dependence of the 
position of the peak of a wave packet. Moreover, some authors have correctly considered 
for a sufficiently complete analysis of the violations of the HE, not only the filter action but 
also the spreading of the wave packets caused by the square-law dependence of the kinetic 
energy on k (sometimes carrying into even negative tunneling times) Jsj and, for some par- 
ticular configurations of two (and more) barriers 45j, the influence of all possible resonances. 
Some additional anomalies with the HE had already been discussed in the last review on 
tunneling time analysis {{J. In spite of quoting the superluminal interpretation present in 
the literature since a long time ago, our discussion concerned with the definition of the 



strict mathematical conditions which limit the ap_ 
for which the Hartman interpretation is valid 



Dlicability of the stationary phase method 



. Strictly speaking, the discussion of su- 
perluminal phenomena is justified only for (ultra)relativistic particles si, Q, for instance 
photons, but not for the tunneling analysis in the case of the non-relativistic Schroedinger 
equation. In fact, up to now, the most interesting experiments concerning the time analysis 
of tunneling processes had been fulfilled with photons. 

In more general lines, there have also been some trying of yielding complex time delays for 
tunneling analysis, ultimately due to a complex propagation constant. In such a framework, 
the supposition of superluminal features is considered artificial since the transmitted peak 
is not causally related to the corresponding incident peak. In certain sense it has caused 
some controversies with denying the physical reality to an imaginary time In parallel 
to the most sensible candidate for tunneling timeshj, J9K a phase-space approach has been 
use to determine a semi-classical traversal time 17|, |4l|, |42j. This semi-classical method 



makes use of complex trajectories which, in its turn, enables the definition of real traversal 



times in the complexified phase space 



41 



421 ]. It is also commonly quoted in the context 



of testing different theories for temporal quantities such as arrival, dwell and delay times 
[7, 9] and the asymptotic behavior at long times 19, 44]. In particular, it suggests that the 
idea of complexifying time should be investigated for some other scattering configurations, 
which reinforces the more general assertion that the investigation of wave propagation across 
a tunnel barrier has always been an intriguing subject which is wide open both from a 
theoretical and an experimental point of view. 
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As a possible solution for partially overcoming some of the incongruities here pointed 
out and quantified, which appear when we compute tunneling phase times in the SPM 
ramework, we have claimed the relevance of the use of the multiple peak decomposition 



391 ] technique previously developed for the above barrier diffusion problem. Essentially, 



we have introduced a way for comprehending the conservation of probabilities for a very 
particular tunneling configuration where the asymmetry presented in the standard case 
was eliminated, and the phase time could be accurately calculated. We mention for a 
subsequent analysis the suggestive possibility of investigating the validity of our approach 
when confronted with the intriguing case of multiple opaque barriers [46], in particular, in 
the case of non-resonant tunneling. Still concerned with subsequent theoretical perspectives, 
the symmetrical colliding configuration also offers the possibility of exploring some problems 
involving soliton structures. To conclude, all the above arguments reinforce the necessity of 
searching the appropriate framework where barrier traversal times can be computed in the 
most generalized way. 
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and the stationary Schroedinger equation, in the presence of a potential barrier 
[48] The validity of the HE was tested for all the other theoretical expressions proposed for the 
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